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Abstract 

We improve, in the energy region between KK tiireshold and ^ 1.4 GeV, 
tiie energy-dependent phase shift analysis of tttt scattering presented in a previous 
paper. For the SO wave we have included more data above KK threshold and we 
have taken into account systematically the elasticity data on the reaction rnr 
KK. We here made a coupled channel fit. For the DO wave we have considered 
information on low energy parameters, and imposed a better fit to the /2 resonance. 
For both waves the expressions we now find are substantially more precise than the 
previous ones. We also provide slightly improved D2 and P waves, including the 
estimated inelasticity for the first, and a more flexible parametrization between 1 
and 1.42 GeV for the second. The accuracy of our amplitudes is now such that 
it requires a refinement of the Regge analysis, for s^^^ > 1.42 GeV, which we 
also carry out. We show that this more realistic input produces tttt scattering 
amplitudes that satisfy better forward dispersion relations, particularly for tt^tt'^ 
scattering. 
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1. Introduction 

In a recent paper by two of usl^l (JRP and FJY), that we will consistently denote by PY05, we have presented 
a set of fits to the data on tttt scattering phase shifts and inelasticities, and we also checked how well forward 
dispersion relations are satisfied by the different tttt scattering phase shift analyses (including our own). 
These various sets differ on the values of the SO phase shifts below KK threshold. We found that some of 
the most frequently used sets of phase shifts fail to satisfy forward dispersion relations and we then presented 
a consistent energy-dependent phase shift analysis of tttt scattering amplitudes that satisfies well forward 
dispersion relations for energies below ~ 1 GeV. Above this energy, we found a certain mismatch between 
the real parts of the scattering amplitudes, calculated from phase shifts and inelasticities, and the result of 
the dispersive evaluations; particularly for 7r°7r° scattering. This we attributed to imperfect experimental 
information in the region IGeV < < 1.4 GeV. 

In the present paper we improve our analysis of the SO wave, the DO wave and, to a lesser extent, 
the D2 and P waves^ in the energy range around and above KK threshold (for the DO wave, we also slightly 
improve the low energy region). For the SO wave we take into account systematically the elasticity data from 
the reaction tttt KK; for the DO wave we include information on low energy parameters, and we improve 
the fit to the /2(127G) resonance, to describe better its width and inelasticity. These two paramctrizations 
are more accurate than what we had in PY05; not only in that they include more data, but also because 
they have smaller errors. 

A slight improvement for the P wave (using a more flexible parametrization) between KK threshold 
and 1.42 GeV is also presented and, for the D2 wave, we improve on PY05 by including its estimated 
inelasticity above ^ 1 GeV. 

Wc have also found convenient to reconsider the Regge analysis, for the energy region above 1.42 GeV, 
particularly in view of the accuracy of the present parametrizations. This we do by taking into account more 
precise values for the intercepts 0:^(0) and ap'{0) than those used in PY05. Although the changes this 
induces are very small, and indeed quite unnoticeable below 1 GeV, the verification of the dispersion relation 
for exchange of isospin 1 above KK threshold is sensitive to this Regge improvement. 

We then also show that, with this more accurate input in the phase shift analysis, the forward tt^tt^ 
dispersion relation is much better satisfied than with the amplitudes in PY05; particularly for energies above 
1 GeV. The Tr^Tr"*" dispersion relation is also improved, but only a little. Finally, the dispersion relation 
for exchange of isospin unity is practically unchanged below 1 GeV, and deteriorates slightly above. The 
new, improved parametrizations, therefore provide a very precise and reliable representation of pion-pion 
amplitudes at all energies: the average fulfillment of the dispersion relations is at the level of 1.05 cr, for 
energies below 0.93 GeV, and of 1.29 ct for energies up to 1.42 GeV. 

2. The SO wave at high energy 

In PY05 we provided fits to data for the SO wave that satisfied forward dispersion relations reasonably well 
below 0.925 GeV, as well as an improved parametrization constrained to satisfy forward dispersion relations 
below this energy and to fit data. Here we will concentrate on a parametrization at higher energies, taking 
care to match it to the low energy one, which we do at 0.92 GeV. 
The information on the SO wave at high energy 

(si/2 > 0.92 GeV) comes from two sources: tttt 
scattering experiments!^"^! and, above KK threshold, also from tttt KK scattering. The second provides 
reliable measurements of the elasticity parameter,^ Vo^H^)'- since there are no isospin 2 waves in tttt —^ KK 
scattering, and the tttt — KK coupling is very weak for P and DO waves, it follows that measurements of the 

^ We will use consistently the self-explanatory notation SO, S2, P, DO, D2, F, . . . for the tttt partial waves. 
^ In the present paper we refer to r; as the elasticity, or elasticity parameter. The inelasticity is -^/l — i]'^. 
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differential cross section for tttt KK give directly the quantity 1 — [?7q X ^ith good accuracy, so long as 
the multipion cross section is small; sec the discussion of this below. 

Below KK threshold we fit data between 0.929 GeV and 0.970 GeV from Hyams et al.,!^] Pro- 
topopescu et al.I^l and from Grayer et al.,W as composed in PYOS:'^! 

4°^ (0.929^ GeV^) = 112.5 ± 13° ; 4"' (0.935^ GeV^) = 109 ± 8° ; 

4°^(0.9522 GeV^) =126±16°; ^^(o.geS^ GeV^) = 134 ± 14° ; (2.1a) 

4°^ (0.970^ GeV^) = 141 ± 18° . 

As explained in PY05, these errors cover the systematic uncertainties, which are large. We also add the 
recent data of Kaminski et al.I^l and we include in the fit the value 

4°^(4m^) = 205±8° (2.1b) 

obtained in the constant K- matrix fit of Hyams et al.,Pl which is compatible with the other data used here. 
Finally, wc include two values that follow from the low energy analysis in PY05 from the global data fit (i.e., 
before imposing forward dispersion relations, to avoid correlations with other waves), 

(5^°^(0.9002 GcV^) = 101.0 ±3.7;°, (5*j"'(0.9202 GcV^) = 102.6 ± 4° . (2.1c) 

To fit the data above KK threshold, we notice that analyses based on tttt scattering experiments 
only determine a combination of phase shift and inelasticity and, indeed, different results are obtained for 
the SO wave in the various analyses. For this wave wc only fit data sets whose inelasticity is compatible with 
what is found in tttt KK scattering, ^''^ in the region 4m|- < s < (1.25 GeV)^. This includes the solution^ 

( ) of Hyams et al.,1^1 the data of Hyams et al.^l (or^ of Grayer et slM) and the data of Kaminski et 

al.PI 

For the elasticity parameter, r]^\ we will improve on the analysis of PY05 by including more data^ 
(especially, tttt — > KK data) and being more realistic in the parametrization. First of all, we remark that the 
modulus squared of the SO amplitude for tttt KK scattering is proportional to |(1 — [??o°'']^), provided the 
two-channel approximation is valid. This is known to be the case experimentally for s^/^ < 1.25 GeV for 
such waves as has been measured, and will very likely be also true for our case (as we verify in Appendix B). 
In this range, the tttt — > KK scattering experiments give the more reliable measurements of the parameter 
li^^ ■ Therefore, in the region s^/^ < 1.25 GeV we fit tttt KK data'^l and, among the tttt tttt data sets, 
only those whose inelasticity is compatible with that from tttt KK below ~ 1.25 GeV. This includes the 

data sets of Hyams et al.'^l (or Grayer et al.'^l); the data from ref. 6, solution ( ); and the data of ref. 5. 

We however do not include in the fits the data of Protopopescu et al.,[^i since they are quite incompatible 
with the TTTT KK information. 

A convenient way to fit phase shift and inelasticity is to use the K-matrix formalism. This has the 
advantage over the method of polynomial fits, used in PY05 (see also Appendix B here), that the relations 
that occur at threshold between cJg"-' and rijj^\ given in Appendix A [Eq. (A. 6)] are automatically fulfilled. 
The method, however, presents the drawback that it is not possible to take into account the existence of 
other channels unless one introduces an excessive number of parameters. This is why we present an alternate 
polynomial fit in Appendix B. Fortunately, the fits given in Appendix B show that the contribution of such 
multiparticle channels is rather small; in fact, within the errors of the two channel fit (this smallness is 

^ ( ) is the preferred solution in the original reference. Unfortunately, this reference only provided numbers 

for the statistical uncertainties. We add to these 5° as estimated systematic error, in agreement with an analysis 
similar to that of PY05. 

The data of Grayer et al. in ref 3, and those of Hyams et al. in ref. 2 come from the same experiment. 
^ In all cases we add an estimated error of 0.04 to data that only give statistical errors. 
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FlGURE lA. K-matrix fit to Sf' (solid line and dark area). Dotted lines: the fit in PY05. 



probably due to the fact that, because of its quantum numbers, the first quasi-two body channel that 
contributes is pp). So we would expect that neglecting those other channels will not produce an excessive 
bias, being anyway covered by our uncertainties. 

To perform the fit, we consider 5^^^ and ri^^^ to be given in terms of the K-matrix elements by the 
expressions (cf. Appendix A) 



ki\k2\<ici'K + kiKii ^ ^ „ 
— , s < 4™! 



l + |fc2|A'22 

1 



2ki[Kii + klK22 det K] 



klKf^ - klKl^ + fc^A:|(det K)^ - f 



{klKl^ + klKl^ + fc2fc2(det K)2 + If - AklklKl^ 



s > Am 



K 



(2.2) 



and 



(0) _ l {l + krk2deiK.f + {k,K^^-k2K22? ^ , 
% {s) - y (1 _ j^^k2 det K)2 + (k,Kn + ^^22)2 ' " 



(2.3) 



- 3 - 



-R. KAMINSKI, J. R. PELAEZ AND F. J. YNDURAIN- 



0.75 




0.25 



TT" 



r 



-- PY05 PolynAmial fit 

T Cohenetal. l-'tt^KK 

• Etkin et al. 

■ Wetzel et al. 

i Hyams et al. 

o Kaminski et al. 
K-matrix up to 1420 MeV 



1000 



1100 



1200 



s^'^MeV) 



1300 



1400 



Figure IB. K-matrix fit to rjj^^ using (2.4), with error given by the shaded area. The dotted lines 
represent the central values and error limits of the old fit in PY05. 



Then we write a standard diadic expansion for K, with a constant background (like, e.g., in Hyams et al.I^l): 

and /i is a mass scale, that we take = 1 GcV. The powers of fi have been arranged so that the a^, (3i, 
jij are dimensionless; they are also assumed to be constant. The pole at simulates the left hand cut 
of K, and the pole at M| is connected with the phase shift crossing 270° around 1.3 GeV; both poles are 
necessary to get a good fit. 

We fit simultaneously all data, above as well as below KK threshold for the phase shift. For rjl^\ 
we also fit all data, tttt —f tttt and tttt —^ KK, over the whole range, which is justified since we are neglecting 
other inelastic channels. We require perfect matching with the lower energy determination of the phase shift 
at 0.932 GeV, as obtained in PY05. We find a x^/d.o.f.=0.6 and the values of the parameters are 

Qi =0.727 ±0.014, a2 = 0.19 ±0.04, /3i = 1.01 ± 0.08; /Ja = 1.29 ± 0.03, 
Ml = 910.5 ± 7 MeV, Ah = 1324 ± 6 MeV; 
711 =2.87 ±0.17, 712 = 1.93 ±0.18, 722 = -6.44 ± 0.17; 

4"^ ((0.92 GeV)2) = 103.55 ± 4.6° . 

Note that Mi indeed lies near the beginning of the left hand cut for KK tttt scattering, located at 
952 MeV. 

The parameters in (2.4b) are strongly correlated. In fact, we have verified that there exists a wide 
set of minima, with very different values of the parameters. This is not surprising, since we do not have 
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sufRciently many observables to determine the three Kij on an energy independent basis. Nevertheless, the 
corresponding values of and ryg"^ vary very little in all these minima, so that (2.4b) can be considered 
a faithful representation of the SO wave for tttt scattering, albeit very likely with somewhat underestimated 
errors due to our neglecting other channels, like tttt — » 47r. The corresponding phase shift and elasticity 
parameter are shown in Figs. lA, IB. Both phase shift and elasticity clearly improve what wc had in PY05. 

It is mainly because of the use of the phase shift in (2.4a) and, above all, the smaller inelasticity 
driven by tttt — > KK data, that we find a substantial improvement in forward n^n^ dispersion relations 
above 1 GeV (see below), as already remarked in PY05. 

3. The improved DO, D2 waves 

3.1. The DO wave 

The experimental data on the DO wave are of poor quality; different experiments give very incompatible 
results, particularly below 0.93 GeV. Above 1.1 GeV, and although the data of refs. 2, 6 (but not ref. 3) are 
compatible, it is better to use directly the Particle Data Table's information!^! on the /2 resonance, which 
summarizes the existing experimental data. 

The reliable information on this wave is then of three kinds. Firstly, in the range around 1.27 GeV, 
we have the referred very precise measurements of the /2(1270) resonance parameters, which give!*' a mass 
M/, = 1275.4±1.2 MeV, a width T/^ = 185.1 ± 3.4 MeV and a tttt branching ratio of 84.7 ± 2.4%. Secondly 
the Froissart-Gribov representation allows an accurate determination of the scattering length, aj"', and 
effective range parameter, 62°^ as shown in PY05;^ one finds'^l 

4°^ = (18.7 ± 0.4) X 10"'* M-^ 6^°^ = (-4.2 ± 0.3) x 10"^ M"^. (3.1) 

This helps us to fix the phase shift at low energy. And thirdly, we have the 1973 data of Hyams et al.Pl, 

Protopopeseu et al.,M and solution ( ) of Hyams et al.'^^ in the range 0.935 GeV < s^^^ < 1.1 GeV, 

which are reasonably compatible among themselves; see Fig. 2. We denote them by, respectively, H73, P 
and H( ). The data we include thus are 



E, in GeV 


cot 


source 


0.935 


10.4 ±2 


P 


0.950 


10.2 ±5.5 


H73 


0.965 


7.3 ±2 


P 


0.970 


4.0 ±2.2 


H73 


0.990 


4.8 ± 1.4 


H73 



(3.2a) 



The values given below in (3.1) are those obtained in PY05, with the old parametrizations. Wc have verified that 
they do not change, within the accuracy of (3.1), if recalculating the Froissart-Gribov representation with the 
parametrizations in the present paper. This of course occurs because the new parametrizations only change the 
amplitudes significantly above 1 GeV, a region to which Oj^^ and 62°^ are almost not sensitive. 
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FlGURE 2. The DO phase shift as determined here (continuous hue; the error is like the thickness 
of the hne) and that from PY05 (broken line) . The experimental data points are also shown. Note 
that the high energy fit is tightly constrained by the /2(1270) mass and width. 



and, above KK threshold, 



in GeV 


cot 


source 


E, in GeV 


cot 


source 


1.00 


5.1 ± 2 


P 


1.05 


3.8 ±0.9 


H73 


1.01 


3.6 ±0.8 


H73 


1.06 


2.8 ±0.8 


H( 


1.02 


3.6 ±0.8 


H( ) 


1.07 


2.5 ±0.4 


H73 


1.03 


3.0 ±0.6 


H73 


1.09 


2.7 ±0.45 


H73 


1.04 


3.3 ±0.8 


P 


1.10 


2.1 ±0.8 


H( 



A few words must be said about the errors in (3.2). Since H( ) do not give errors, we take them as 

equal to those of P. We also multiplied all errors by a factor 2, to take into account the estimated systematic 
errors (for e.g. P, estimated as the difference between the fits XIII and VI in ref. 3). 

We present the details of the fits. To take into account the analyticity structure, we fit with 
different expressions for energies below and above KK threshold, requiring however exact matching at 
s = Am\. Below KK threshold we take into account the existence of nonnegligible inelasticity above 
1.05 GeV, which is near the wtt or pirir thresholds, by choosing a conformal variable w appropriate to a plane 
cut for s > (1.05 GeV)^. So we write 



1/2 

(3.3a) 

= 1.05 GeV. 



cot 4"' (s) = 1^ (M|^ - s) Af 2 + Biu;} , s< ^m" 
fs — vs 
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— DO fit 

O fi-omf,(1275) inPDG 

* Hyams et al. (— ) 

■ Piotopopescu et al. (Table XIII) 

• Protopopescu et al. (Table VI) 
□ Hyams et al ("73) 
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Figure 3. Fit to r?^°^ (continuous line 
and dark area that covers the uncertainty), 
from PY05. Data from refs. 2, 3, 4, 6. 
The elasticity on the /2(1270), from the 
PDT,'*' is also shown (large white dot). 



The mass of the /2 we fix at Mf^ — 1275.4 MeV; no error is taken for this quantity, since it is negligibly 
small (1.2 MeV) when compared with the other errors. We fit the values of 0^2^ and 62*''' given in (3.1) and 
the data in (3.2a). We find the values of the parameters 

Bo = 12.47 ±0.12; Bi = 10.12 ± 0.16. (3.3b) 

We note that the series shows good convergence. 

Above KK threshold we use the following formula for the phase shift: 

cot 4°^ (s) = 1^ (M/, - s) Ml \ Bm + Bhiw\, s> Am\ ; 

r / (3.4a) 

^^ 7i-v^ . ,1/2^ 1.45 GeV. 

This neglects inelasticity below 1.45 GeV, which is approximately the pp threshold; inelasticity will be added 
by hand, see below. We then fit the values for the width of /2(1270), as given above, and the set of data in 
(3.2b). We get 

Bm = 18.77 ±0.16; B/a = 43.7 ± 1.8. (3.4b) 
As stated, we have required exact matching of high and low energy at KK threshold, where our fits give 
cot Jj"'' (^tit-'k) = 4.42 ± 0.04. This matching implies that there is a relation among the four BiS, so there are 
in effect only three free parameters. The overall chi-squared of the fit is very good, /d.o.f. = 9.8/ (18 — 3). 

We note that, although not included in the fit, our new DO phase shift fits better than the old PY05 
one the data points of Hyams et al.,'^' solution ( ), above the /2 resonance: see Fig. 2. 

The data for the inelasticity are not sufficiently good to improve significantly the fit in PY05; so we 
simply write, as in ref. 1, 

r 1, s< 4m|., 

(^)= <^ l_e-Mf) , s>Aml; 6 = 0.262 ± 0.030, (3.4c) 
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This probably only provides a fit to elasticity parameter on the average, but we have not been able to find 
a clear improvement on this. The corresponding elasticity parameter is shown'' in Fig. 3. 

Two important properties of the new fit are that it reproduces better than the one in PY05 the 
width and inelasticity of the /2 resonance, which is the more salient feature of the DO wave, and that it 
is more precise than what we had in PY05. This improvement of the DO wave, although it docs not give 
a phase shift very different from that in PY05, also contributes a nonnegligible amount to the improved 
fulfillment of the w^w^ dispersion relations. 

3.2. The D2 wave 



In PY05 we fitted the D2 wave with a single parametrization over the whole energy range up to 1.42 GeV, 
and neglected inelasticity. We wrote 

cot Si'\s) = ^ {So + B,w{s) + BMsf} ^^Mffl^) - s ^^"^"^ 
with A a free parameter fixing the zero of the phase shift near threshold, and 

Since the data on this wave are not accurate we included extra information. To be precise, we incorporated 
in the fit the value of the scattering length that follows from the Froissart-Gribov representation (PY05), 

= (2.78 ±0.37) X 10-''M-^ 

but not that of the effective range parameter, 

^(2) ^ (-3.89 ±0.28) X 10"* M^-^. 

We got a mediocre fit, /d.o.f. = 71/(25 — 3), and the values of the parameters were 

Bo = (2.4±0.3) X 10^ Bi = (7.8 ± 0.8) X 10^ B2 = (23.7 ± 3.8) x 10^ Z\ = 196 ± 20 MeV . (3.5b) 

(2) 

The corresponding numbers for the scattering length and for the effective range parameter 62 that follow 
from this are 

af^ = (2.5 ± 0.9) X 10"* M^'^; b^^^ = (-2.7 ± 0.8) x 10"* M^"^. 

The last is a bit away from what one has from the Froissart-Gribov representation, but still is compatible 
at the 2 a level. The low quality of the fit may be traced to the fact that the various data sets are not very 
compatible among themselves. Therefore, there is no chance to improve the fit as we did for the DO wave 
(where wc; had the very precise data on the /2 resonance). We here merely improve the treatment of this 
wave by including the inelasticity by hand. 

To get an estimate of the inelasticity we have two possible methods: we can take the inelasticity 
to be similar to that of the DO wave; or we can make a model calculation. For example, that of ref. 9, 
in which the authors assume inelasticity to go via rho intermediate states, fixing the coupling parameters 
to reproduce the properties of the better known waves. Both methods yield negligible inelasticity below 
pTTTT threshold, and something around 5% inelasticity at the highest energy considered, 1.42 GeV. For the 
elasticity parameter we thus simply write, above 1.05 GeV, 

ri!^\s) ^ 1- e{l- s/sf, s^/^ ^ 1.05 GeV, e = 0.2 ± 0.2; (3.5c) 

this is negligible up to 1.25 GeV and, above that, covers both what was estimated in ref. 9, and the fact that 
experiments fail to detect inelasticity. 



Although there is nothing new in this fit, we show the picture because we had not shown it in PY05. 
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FlGURE 4. The fit to the P wave above KK threshold (continuous line and dark area that covers the 

uncertainty), with data from solution ( ) of Hyams et alJ^' and of Protopopescu et al.'^' Note that 

the errors shown for the data are only the statistical errors; systematic errors, estimated as in the text, 
about double them. The broken lines are the phase shift and elasticity parameter of PY05. The effect 
of the <j!>(1020) resonance is not shown in this figure. 



4. The improved P wave between KK threshold and 1.42 GeV 



We next fit the P wave above 2mK — 0.992 GeV, incorporating in the fit the data from solution ( ) of 

Hyams et al.,'^' besides the data from Protopopescu et al.W (the last is the one more compatible with what 
one finds from the pion form factor). We have added estimated errors of 2° to the phase shift and 0.04 to 

the elasticity parameter for the data of solution ( ) in Hyams et al.,'^' since no errors are provided in 

this reference. We now use one more parameter both for the phase shift and for the elasticity parameter 
than what wc had in PY05, writing 

5i(s) = Ao + XiiJs/Amj, - 1) + X2{Js/4mj, - 1)^, 

^ ^ (4.1) 



r]i{s) — l~ei^l — imj^/s — £2(1 — 4rn|-/s); s > 4m 



The phase at the low energy edge, (5i (0.992^ GeV^) = 153.5 ± 0.6° , is obtained from the fit to the form 
factor of the pion (ref. 9; see also ref. 1). This fixes the vale of Aq. 

The fits are reasonable; we get /d.o.f. = 0.6 for the phase and /d.o.f. = 1.1 for the elasticity. 
We find the parameters 

Ao =2.687 ±0.008, Ai = 1.57 ± 0.18, A2 = -1.96 ± 0.49; , , 

(4.2) 

ei =0.10±0.06, £2 = 0.11±0.11. 
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The only noticeable differences with the fit in PY05 is that the phase shift and elasticity parameter are now 

less rigid, that wc match the low and high energy expressions at KK threshold, and that the inelasticity is 
now somewhat larger than what we had in PY05. This improved solution, together with that in PY05, are 
shown in Fig. 4. 

Another matter is the contribution of the (^(1020) resonance. This can be included in the standard 
way, by adding to the P wave a resonant piece 



where /i(s) is normalized so that, in the elastic case. 

/i = sin,5ie'^^ 

and 



(4.3a) 





k2 


3 


sV2 


_fc2(M|)_ 





-Btt 



s > 4m 



.sV2 



(4.3b) 



here k2 = a/s/4 — m|-, and the width and tttt branching ratio of the 0(1020) resonance are = 4.26 ± 
0.05 MeV and i?7r7r = (7.3 ± 1.3) x 10^^. Something similar could be done for the contribution of the oj. 

The influence of these resonances is totally negligible and, in fact, we will not include them in our 
calculations of dispersion relations below. 



5. Improvement of the Regge input 



To evaluate the dispersion relations we need an estimate for the high energy (s^/^ > 1.42 GeV) scattering 
amplitudes. This is furnished by the Regge model. We have here three amplitudes, one for each of the 
exchange of isospin 0, 1 and 2. Wc first take, for these Regge amplitudes, the results of the fits in ref. 11; 
see also PY05, Appendix B. Then we will consider improvement of the Regge parameters. 

The expressions for the amplitudes for exchange of isospin 1 and are, respectively, 

IinFUt=i)(s,0) ~ /3p(0)(s/s)"''W, s>(1.42GeV)^ (5.1a) 



and 

ImF^^*=°)(s,0) ~ P{s,0) + P'{s,0), s > (1.42 GeV)2; 



(5.1b) 



P(s,0) =/3p(s/s)«^(°), P'{s,0) = pp, (s/s)«^'W. 
In both expressions s = 1 GeV^. The values of the parameters are (ref. 11 and ref. 1, Appendix B) 

/?p(0) = 1.02 ±0.11; ap(0) = 0.52 ±0.02 (5.2) 

and 

/5p = 2.54 ±0.03, /3p' = 1.05 ±0.02; ap. (0) = ^^(0), ap(G) = 1. (5.3) 
For exchange of isospin 2, which is very small, we also take the amplitude of ref. 11: 

ImF(^*=2)^_g^0) ~ /32(s/s)^"'■^°^"^ /32 = 0.2±0.2; s > (1.42 GeV)^ (5.4) 

The first two, however, will now be improved: as we have seen in previous sections, the precision of 
our new parametrizations in the intermediate energy range (~ 1 to ^ 1.4 GeV) is such that one is sensitive 
to small details of the Regge amplitudes; so, it is convenient to re-assess the derivation of the values for the 
Regge parameters in Eqs. (5.2) to (5.3). 
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The expressions (5.2), (5.3) were obtained in ref. 11 and PY05 as follows. We fixed 0:^(0) as the 
average between what is found in deep inelastic scattering, [^^1 ap = 0.48, and in the analysis of hadron 
collisions by Rarita et al.,'^'^1 who get ap = 0.56. We also imposed degeneracy, so that the intercept of p and 
P' were forced to be the same. We then fitted experimental tttt cross sections, which gives /3p(0) = 1.0 ± 0.3, 
and improved this result demanding fulfillment of a crossing sum rule. For isospin zero exchange, the 
expression (5.3) was obtained requiring simultaneous fits to tttt, -kN and NN data, using factorization, 
fixing the intercept of the P' to 0.52 (as already stated). 

However, more complete fits''^'*! than that of Rarita et have been performed in the last years; 

especially, for the rho trajectory, individual data on pp, pp and np have been included in the fits, which 
permits improvement of the determination of the rho parameters using factorization. These fits, in particular, 
allowed a relaxation of the exact degeneracy condition 0:^(0) = ap'(O), and yield central values for the rho 
intercept ap = 0.46, more in agreement with the result from deep inelastic scattering. For ap' one finds a 
value higher than for the rho intercept: api = 0.54. 

We may then repeat the analysis of ref. 11, but fixing now the intercepts of rho and P' trajectories 
to the likely more precise values 

ap(0) = 0.46 ±0.02, ap/ (0) = 0.54 ± 0.02, (5.5a) 

with conservative errors. We also here improve the error estimate for the rho residue (}p with the crossing sum 
rule, as we did in PY05 to get (5.2). This sum rule we calculate using the new phase shifts and inelasticities 
we have evaluated in the present paper. We then find, 

Pp = 1.22 ± 0.14, /Jp = 2.54 ± 0.04, /3p/ = 0.83 ± 0.05. (5.5b) 

The errors arc slightly larger now, which is due to the fact that we do not impose the exact degeneracy 
relation ap(0) = ap'(O). For the amplitude with exchange of isospin 2, we still keep (5.4) since no new 
information is available. 

The difference between what we have now, (5.5), and what was used in PY05 is much smaller than 
what would appear at first sight; in fact, because the a(0) and /3 are strongly correlated, the changes in one 
quantity are compensated by those in the other: the amplitudes described by (5.5) and (5.2, 3) are very 
similar in the energy region of interest (cf. Fig. 5). However, these amplitudes differ in some details. So, 
the rho amplitude described by (5.5) is tilted with respect to that given by (5.2): the amplitude described 
by (5.5) is slightly larger than that described by (5.3) below 5 GeV, where they cross over, and is larger 
above this energy. Likewise, for exchange of isospin zero (5.5) gives a smaller amplitude at low energy, which 
then crosses over the amplitude given by (5.3) at higher energy. 

As just stated, these changes induced by using (5.5) do almost compensate each other and, indeed, 
they have only a minute effect in dispersion relations below 1 GeV. At the level of precision attained by our 
parametrizations in the region above 1 GeV, however, the dispersion relations are sensitive to the details of 
the Regge behaviour; because of this, we will evaluate the dispersion relations with both (5.2, 4) and with 
(5.5). 

A last question related to the high energy, s^/-^ > 1.42 GeV, input is the matching of the Regge 
amplitudes to the amplitude obtained below 1.42 GeV with our phase shift analyses. Although we have 
verified that the low and high energy amplitudes are compatible, within errors, at 1.42 GeV, we have not 
required exact matching. The reason for this is that at the lower energy Regge range, say below ~ 1.8 GeV, 
some amplitudes still present structure; for example, for the amplitudes with isospin unity, the structure 
associated with the p(1450), /9(1700) and /33(1690) resonances. It is true that these resonances couple weakly 
to TTTT, but, at the level of precision required in the present paper, this is not negligible: as happens in the 
case of TT+j) scattering (see e.g. Fig. 2 in ref. 11), one expects the Regge amplitude to provide only a fit in 
the mean. This mismatch produces distortions near the boundary, s^/^ = 1.42 GeV, clearly seen in some 
of the dispersion relation calculations below; particularly, for tt^tt''" scattering, where the P and, to a lesser 
extent, the F waves are important. We have done nothing to correct this distortion which, anyway, only 
affects the points very near 1.42 GeV. The alternate possibility, which would be to use phase shift analyses 
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Figure 5A. 

The scattering amplitude 
Im F'^'^^-* (s, 0) as described 
by (5.2), broken line, and (5.5), 
solid line with error included 
(gray band). 
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Figure 5B. 

The scattering amplitude 
ImF(^*=°)(s,0) as described 
by (5.3), broken line, and (5.5), 
solid line with error included 
(gray band). 



up to higher energies, say 1.8 GeV, would only make matters worse since it would have to contend with the 
nonuniqueness and unreliability of the experimental data in that region, as discussed for example in rcf. 15. 
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6. Forward dispersion relations 



In this Section wc will evaluate forward dispersion relations for the three independent tttt scattering am- 
plitudes. For these calculations we will take the parameters for all partial waves from the fits to data* in 
ref. 1 (PY05), except for the SO and P waves above 0.92 GeV, where we use the expressions found in the 
present paper (for the SO wave, with the K-matrix fit), and for the D2 wave, where wc take into accoimt the 
inelasticity above 1.05 GeV. For the DO wave we use the expressions given in the present paper all the way 
from threshold. 

To measure the fulfillment of the dispersion relations we calculate the average chi-squared, . This 
is defined as the sum of the squares of the real part minus the result of the dispersive integral, divided by the 
(correlated) errors squared; this we do at energy intervals of 25 MeV, and divide by the number of points. 
Note however, that this average docs not come from a fit to the dispersion relations, but is simply a 
measure of how well the forward dispersion relations are satisfied by the data fits, which are independent for 
each wave, and independent of dispersion relations. When calculating this y^, we first use the parameters 
for phase shifts and inelasticities in PY05; then, wc replace the relevant waves by the ones in the present 
paper; and, finally, we also replace the PY05 Regge parameters with the ones in Eq. (5.5). 

6.1. The TT^TT*^ and tt'^tt"'" dispersion relations 



We first evaluate the forward dispersion relation for 7r°7r° scattering, the one that was worse verified in PY05 
and the one for which the improvement due to the new parametrizations is more marked. We write 

The result of the calculation is shown in Fig. 6, where the continuous curve is the real part evaluated from 
the parametrizations, and the broken curve is the result of the dispersive integral, i.e., the right hand side 
of (6.1). 

The fulfillment of this dispersion relation improves substantially what we had in PY05:^ the changes 

in the average chi-squared are 

7r°7r° : PY05 New phase sh. New Regge 

^2 = 3.8 ^ 1.52 ^ 1.41, for s^/^ < 930 MeV, (6.2) 

^2 = 4,8 ^ ^ 163, for s^/^ < 1420 MeV. 

Here and in similar expressions below, "New phase sh." means that we use the new, improved phase shifts 
(and inelasticities) of the present paper; "New Regge" means that we also use the new Regge parameters in 
(5.5). In both cases we use the K-matrix fit for the SO wave, Eqs. (2.4). 

The improvement obtained for tt'^tt" when using the new phase shifts is more impressive if we 
remember that the errors we have now for the SO wave above 0.92 GeV, and for the DO wave in the whole 
range, are substantially smaller than what we had in PY05. It is also noteworthy that the improvement 
in the dispersion relation is due almost exclusively to the use of the new phase shifts and inelasticities in 



In PY05 wc gave two sets of phase shifts and inelasticities: one by fitting directly the various sets of experimental 
data (Sect. 2 in ref. 1); and a set obtained by requiring, besides fit to data, fulfillment of dispersion relations 
(summarized in Appendix 1 of ref. 1). In the present paper we of course only use the amplitudes obtained in PY05 
by fitting data, since the ones improved with dispersion relations use a high energy {s^^^ > 0.92 GeV) input that 
is superseded by our calculations in the present paper. 

Of course, we here compare with the results obtained using the fits to data, before improving them by requiring 
fulfillment of the dispersion relations at low energy 
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Figure 6A. 

The tt"tt^ dispersion relation with 
the new SO, DO and D2 waves. 
Continuous line: real part, evalu- 
ated directly with the parametriza- 
tions (the gray band covers the er- 
ror). 

Dashed line: the result of the dis- 
persive integral, with the Regge pa- 
rameters of (5.5). 
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Figure 6B. 

The tt^tt" dispersion relation with 
the old, PY05 SO, DO and D2 waves. 
Continuous line: real part, evalu- 
ated directly with the parametriza- 
tions. 

Dashed line: the result of the dis- 
persive integral, with the Regge pa- 
rameters as in PY05. 



the range ^ 1 to 1.42 GeV; the improvement due to introducing the Regge behaviour (5.5) is much more 
modest. 

The dispersion relation for tt'^tt^ scattering reads, with Fq^{s) the forward tt'^tt^ amplitude, 

{2s' -AM^)lmFo+{s') 



ReFo+(s) - Fo+(UO = ^P.P. / ds' - , 

^ TT s'(s'-.s)(s'-4M2)(s' + .s-4M2) 



(6.3) 



In Fig. 7 we show the fulfillment of (6.3), both with what we had in PY05 and with the new phase shifts 
and Regge parameters. 
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Figure 7A. 

The Tr^Tr"*" dispersion relation with the 

new P and D2 waves. 

Continuous line: real part, evaluated 

directly with the parametrizations. 

Dashed line: the result of the dispersive 

integral, with the Regge parameters of 

(5.5). 
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Figure 7B. 

The TT^Tr"*" dispersion relation with the 
old, PY05 P and D2 waves. 
Continuous line: real part, evaluated 
with the parametrizations. 
Dashed line: the result of the disper- 
sive integral, with the old, PY05 Regge 
parameters. 



The forward dispersion relation for 7r°7r+ scattering was already very well satisfied with the param- 
eters in PY05; it becomes slightly better satisfied now. The changes in the average chi-squared are 



X 



PY05 
= 1.7 
= 1.7 



New phase sh. 
1.75 
1.60 



New 
^ 1.60, 
1.44, 



for 5^/2 < 930 MeV, 
for 5^/2 < i_42o MeV; 



(6.4) 



The improvement here, although existing, is rather small: not surprisingly as the corresponding 
amplitude does not contain the SO or DO waves. The amelioration is due only to use of the new Regge 
parameters from Eq. (5.5). 

The fact that both the dispersion relations for tt'^tt" and tjP'k^ improve with the present parameters 
for the P' trajectory confirms the correctness of the procedure for determining it which we developed in 
Sect. 5. 
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Figure 8 A. The dispersion relation 
for the It — 1 amphtude, calculated 
with the new amplitudes. 
Continuous line: real part and error 
(shaded area) evaluated directly with 
the parametrizations. 
Dashed line: the result of the disper- 
sive integral, with the Regge parame- 
ters given in (5.5). 




PY direct from data 
PY dispersive from data 



400 



600 



800 
s''^ (MeV) 



1000 1200 1400 



Figure 8B . The dispersion relation 
for the = 1 amplitude with the old, 
PY05 SO, DO and P waves. 
Continuous line: real part and error 
(shaded area) evaluated directly. 
Dashed line: the result of the dispersive 
integral, with the PY05 Regge param- 
eters for the rho. 



6.2. The dispersion relation for the It = 1 scattering amplitude 



The dispersion relation for the It — ^ scattering amplitude does not require subtractions, and reads 

2s~Ull „„ . , IniF(^*=i)(s',0) 



ReF('^'=i)(s,0) = 



P.P. 



4M? {s' - s){s' + s - AMI)- 



(6.5) 



The result of the calculation is shown in Fig. 8. 

In this case the contribution of the Regge piece is very important, although the details only matter 
in the region above 1 GeV. Here the fulfillment of the dispersion relation becomes entangled with which 
Regge behaviour one uses; particularly since we now have SO and DO amplitudes with very small errors 
above 1 GeV, which is where the detailed shape of the Regge amplitude has more influence. 
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The changes in the from what we had in PY05 are 
It = 1 : PY05 New phase sh. New Regge 
= 0.2 ^ 0.57 ^ 0.32 

= 1.4 2.32 1.76 



for s^/^ < 930 MeV, 
for s^/^ < 1420 MeV . 



(6.6) 



The conventions are Hke in (6.2) above. 

The dispersion relation deteriorates a little, which indicates that the rho Regge parameters may still 
be improved. In fact, it is remarkable that the simple change of (5.5) in place of (5.2, 3) improves so clearly the 
dispersion relation above 0.9 GeV for exchange of isospin unity, while leaving it almost unchanged below this 
energy for all processes. This confirms that the Regge parameters are much better determined for exchange 
of isospin zero than for exchange of isospin 1, and indicates that a complete treatment of dispersion relations 
(in particular, using them to improve the scattering amplitudes) may require simultaneous consideration of 
the Regge parameters and of the parameters of the phase shift analyses, as in fact was done in PY05. We 
will leave this for a forthcoming paper, where we will also study the improvement of our parametrizations 
using the dispersion relations as well as Roy equations. 

Finally, we mention here that, although the improvements in the present paper only affected the 
various waves above ^ 1 GeV (with the exception of the very small change of the DO wave below KK 
threshold), there is a systematic improvement of the dispersion relations also below that energy, which is a 
nontrivial test of the consistency of the parametrizations below and above KK threshold. 

7. A brief discussion 

The results of the present article show that, if we improve the scattering amplitudes above ^ 1 GeV us- 
ing more reliable data sets that those we had in PY05, the ensuing amplitudes verify much better forward 
dispersion relations, especially above KK threshold; but also below it. Forward dispersion relations, par- 
ticularly for and TT^Tr"*" scattering, which (as discussed in PY05) have important positivity properties, 
constitute a very stringent filter when used to discriminate against spurious parametrizations or calculations, 
as discussed in PY05 and ref. 15. The fact that, with the small errors we have now, all values for the 
are below the 1.8 level, implies that a small change in the parameters would ensure complete fulfillment 
(within errors). However, it is clear that, although small, some alterations are to be expected of the various 
parameters if we require the amplitudes to verify dispersion relations at the = 1 level, which we will do 
in a forthcoming article. 

These changes are forced by the fact that the dispersion relations are not yet perfectly satisfied. 
With respect to this, we have three suspects here. First of all, we have that the experimental data for the D2 
wave (which contributes to all processes) are of such a kind that our fit cannot be very reliable for the phase 
shift above 1 GeV, and is almost pure guesswork for the inelasticity. In fact, already in PY05 we discovered 
that requiring fulfillment of the dispersion relations, within errors, forces a change by more than 1 cr in the 
phase shift parameters for this D2 wave. 

The second possible culprit is the inelasticity for the DO wave. Although it fits (by construction) 
that of the /2 resonance, the expression we have used is, probably, too rigid. There is unfortunately very 
little one can do here, since the quality of the data does not allow an accurate treatment. 

The final possible culprit is the isospin 1 Regge amplitude: there is perhaps room for improvement 
here. The same is true, albeit to a lesser extent, for the amplitudes for P' and for exchange of isospin 2. 
(Alternatively, it may turn out that, once the D2 wave is improved, any change in the Regge parameters is 
unnecessary) . 

Finally, it is clear that one cannot improve our amplitudes much, since they are quite good to begin 
with. However, and based on the preliminary results that we have at present, we expect to show, in a 



-17- 



-R. KAMINSKI, J. R. PELAEZ AND F. J. YNDURAIN- 



forthcoming article, that it is still possible to hone our amplitude analysis by requiring fulfillment of the Roy 
equations and, especially, of forward dispersion relations over the whole energy range. 

Appendix A. The K-matrix formalism 



The phase shift and elasticity parameter r/ for the SO partial wave, for tttt scattering, are defined as 



/ii(s) =sin5^e'*'', 
r;e^'^- -1 



s < 4m|f ; 



s > Am 



(A.1) 



K- 



Wc have changed a little the notation with respect to the main text; thus, 5^^ is what we called 5^^^ before, 

/ii was called /q°'' in the main text, etc. The index (11) in /n is a channel index; see below. Also, we do 
not write angular momentum or isospin indices explicitly. We assume here, as in the main text, that there 
are only two channels open (which is likely a good approximation below ~ 1.25 GeV, and not too bad up to 
1.42 GeV): 

(11) : TTTT ^ tttt; (12) : tttt ^ KK; (22) : KK KK. 

Because of time reversal invariance, the channels tttt — > KK and KK — > tttt are represented by the same 
amplitude. We then form a matrix, with elements fij, i, j = 1,2, /n = fTrv^v-K, etc.: 



TTTT— tTTTT J-^TT^KK \ _ 

KK IkK^KK 



r]e 



1 



2i 



2i 



(A.2) 



f ={k-l/2K-lk-l/2_i} \ k= 



Sk is the phase shift for KK KK scattering. Below KK threshold, the elasticity parameter is ??(s) = 1; 
above KK threshold one has the bounds < ?7 < 1. 
We write f as 

)• (^•^) 

ki are the momenta, ki = a/ s/4 — M^, k2 = i/s/4 — m\. Then, analyticity and unitarity imply that 
K is analytic in s through the KK threshold; hence, it only depends on A:|: Kij = Kij{k2). This is the 
well-known K-matrix formalism, which the reader may find developed in detail in the standard textbook of 
Pilkuhnl^^l or, perhaps more accessible, in the lecture notes by one of us;[^^l and, applied to the SO wave in 
TTTT scattering, in ref. 2. 

Because of (A.2, 3), one can express and 7] in terms of the K^ as 

' ki\k2\deiY. + kiKii 



tan^TT 



1 + \k2\K22 
1 



s < 4m|-, 



2ki[Kii + klK22deiK\ 



klK'h - klKl^ + klkl{det Kf - 1 



(A.4a) 



{klKl^ + klKl^ + kfk^idet K)2 + 1)2 _ Akfk^Kf^ }, s>4m 



if 1 



one also has 



/(l + fcifc2detK)2 + (/eii^n-fei^22)" 2 



(1 - kik2 det K)2 + {kiKn + k2K22Y ' 
The sign in the surd in (A. 4a) is to be taken positive if, as happens in our case, Kii{k2 = 0) > 0. 



(A.4b) 
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From the relation between the phase shift above and below threshold, and also with the elasticity, it 
may appear that one could write an expansion for 5.,r{s) below threshold and from it, deduce corresponding 
expressions for 777r(s) and for 57r(s) above threshold. This conies about as follows. Let us define 5^{s) to be 
the phase shift below threshold, and 5%{s) that above threshold, both as given in (A. 4a). Write the Taylor 
expansion 

oo 

5l{s) = ^ a^K^'lm^, Oo = do (A.5) 



and K = \k2\ below threshold. Substituting this into the expression, valid below threshold, 

and continuing this across the cut in the variable k = —\k2 above the threshold we find the expression, valid 

for s > 4:m'^, 

Q2{aik2/mK-a3kl/m% + ---)Q2i(do-a2kl/mj^+a4k2/m'^ + ---) _ ^ 

On comparing with the expression above threshold given in (A.l) we find 

S^{s)=dQ-a2kl/m% + aik^/m%-\ ; 

,7(s) =e2(«i'==/™^-«3'=2Vml.+-)^ {s>4mj^). 

However, the convergence of (A. 6) can only be guaranteed in a disk touching the left hand cut of the K- 
matrix, a cut due to the left hand cut in KK nn scattering, '^'^ that runs up to s = 4(m^ — M^): therefore, 
only for |A;2| < M^. From a practical point of view, we have checked numerically that fitting with (A.5), (A. 6) 
represent reasonably well S^{s) and S^{s) (this one with irrealistic errors) but does certainly not represent 
T){s), in the region away from s = 4?77,|-, unless one adds an inordinately large number of parameters. 

On the other hand, it is clear that all three S^{s) , S^{s) and 77(5) are continuous functions of, 
respectively, k, fc^ and ^2- Therefore they can be approximated by polynomials in these variables over the 
whole range, even if they are not one the continuation of the other. 

Appendix B. Polynomial fit 



We present here a polynomial fit to phase shift and elasticity parameter in which the three quantities: phase 

shift below KK threshold, phase shift above this threshold, and elasticity are fitted separately. Although 
this fit is less reliable than the K-matrix one, especially near KK threshold, it will allow us to test the 
importance of multibody channels, not taken into account in the K-matrix fit. 
For the phase we write 



4°)(. 



do + a 1- — 5-, 

niK 

k^ 



(0.92 GeVy <s< 4m|.; 



do + B-^, Amji <s< (1.42 GeV)^; k2 = ^/Tfl- 



(B.l) 



m 



m 



K- 



K 



^^\\ is not difficult to check that, although /ii(s) or jii{s) have no left hand cut above s = 0, 7)(s) and <5*(s) do. For 
e.g. the first, we use (A.l) and find 



2 

f] = 



(2i/il + 1) (21/22 + 1) 



(2i/ll+l)(2i/22 + l)+4/22 

from which it is obvious that 77 inherits the left hand cut of /12. 
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FlGURE 9A. Comparison of the fits to : polynomial [Eq. (B.2)] given by the dashed lines, and 
with the K-matrix [Eq. (2.4)] (solid line and dark area). 



The parameters do, a and b are strongly correlated. One can get parameters with low correlation by elim- 
inating the parameter b in favour of the phase shift di at a low energy point, that we conveniently take 
s^/^ — 0.92 GeV: note that, unlike for the K-matrix fit, we now match low and intermediate energy fits at 
0.920 GeV. We thus rewrite the parametrization as 

k2\ , |fc2p r, , |fc2(0.922GeV2)| 



(0), s I rriK 1/02(0.922 GeV^)|H t^^k 



do + + ,, ,^1;' ,,2m9 {di-do- a'-^ ^ \ , (0.92 GcV)^ < s < 4mj,; 



k. 



2 

do + B^, 4to| <s< (1.42 GeV)2; fcz ^s/4- 



m 



m 



K- 



K 



(B.2a) 

In the previous Appendix A we presented a discussion about these expansions. From it it follows that, while 
the expansion below threshold can be considered as convergent in the range of interest here, 0.92 GeV < 
gi/2 £ 2mK, the expansion above threshold (both for (5g°^ and r]^\ see below) should be taken as purely 
phenomenological. In particular, we do not impose the equality b — ~B that would follow if we took (B.2a) 
to be a Taylor expansion (see Eq. (A. 6) in the Appendix). It is possible to fit requiring b = —B, at the cost of 
adding an extra parameter in (B.l), c\k2\'^/rn^. The resulting fit is not satisfactory: it presents excessively 
small errors for s > ^rri]^, due to the forced relation b = —B, which should only be effective near threshold, 
the only region where the expansion converges. 
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We fit separately data above and below KK threshold. The fit returns a / d.o.f . = 0.4 below 
threshold, and x^/d.o.f. = 0.9 above threshold; the values of the parameters are 

rfo = 218.3 ± 4.5° , a = -537 ±41°, rfi = 4°^(0.9202 GeV^) = 102.6 ± 4° (B.2b) 

and 

B = 96 ± 3° . (B.2c) 

The resulting phase shift is shown in Fig. 9A, compared with the K-matrix fit. 

Above 1.25 GeV, the two channel formalism is spoiled by the appearance of new channels, notably 
TTTT 47r, so one does not have an exact connection between the data on tttt KK and riQ \ In fact, 
the numbers one gets for t^q"' from tttt — > tttt, and those that follow from tttt — > KK , assuming only two 
channels, are slightly different; see below. We may take this into account by using a polynomial fit (instead 
of a K-matrix one, as we did in the main text). 

We next make a polynomial fit to the elasticity parameter writing 

''r=l-(^.^ + 4 + «i)^ (B.3a) 

In principle, the values of the arc related to the a, c, . . . of (B.l); see Appendix A, Eq. (A. 6). However, we 
will not impose such relations, but will consider the as phcnomcnological parameters, completely free. The 
reason is that the expansion 2(ak2/'m,k+ck2/m^ + - • •) is very poorly convergent above ^1.2 GcV: something 

that is a disaster for rjl^^ , since the expansion appears in an exponent (the reason for this divergence, that 
can be traced to the left hand cut in tttt KK scattering, may be found in Appendix A). Therefore, we 
would need to add extra phcnomcnological terms, very large, to compensate for that: it is more reasonable to 
make the fit phcnomcnological from the beginning. What we lose by so doing is that we are overestimating 
the value of r]j^^ (s) for s very near threshold, say for 2mK < s^^^ ^ 0.997 GeV, a reasonable price to pay 
to get a good description of the elasticity in the rest of the range. 

We will consider the following possibilities: (a) To fit only tttt nn data above 1.25 GeV, in 
principle the more reliable option; (b) To fit aiso nn — > KK data above 1.25 GeV; (c) To fit only tttt — > KK 
data, in the whole range. Of course, below 1.25 GeV we include both tttt — > KK and tttt — > tttt data in the 
fits (a) and (b). We find 

ei =5.45 ±0.04, 63 = -30.0 ± 0.15, £3 = 46.3 ± 0.5; xVd.o.f. = 1.1 (a); 
ei =5.27 ±0.04, ea = -28.2 ± 0.15, 63 =42.2 ±0.5; xVd.o.f. = 1.1 (b); 
ei =5.77 ±0.05, e2 = -32.9 ± 0.2, £3 = 51.1 ±0.5; xVd.o.f. = 0.2 (c). 

Note that the errors given here are purely nominal, as the parameters are very strongly correlated, while 
they were here treated as uncorrelated. Note also that the three fits arc less separated that it would seem, 
precisely because of that correlation. Finally, we remark that the value of rj^y'^ that follows from (c) is larger 
than what follows from (a) or (b); and (b) also slightly above (a). These two features constitute very nice 
consistency tests, since taking rjl^^ to be given from tttt KK as if only two channels were present must 
surely underestimate the inelasticity; particularly above ~ 1.2 GeV, where the process tttt — > 47r is expected 
to become nonegligible. 

We have verified that one may cover the; two fits (a) and (b) (and even overlap (c), at the edge of 
the error region) by taking as central value that of the fit (a) above and slightly enlarging the errors. We 
then get our best result: 

ei = 5.45 ± 0.06, e2 = -30.0 ± 0.2, £3 = 46.3 ± 0.8; (B.3b) 

the errors may now be taken as uncorrelated. The resulting elasticity may be seen in Fig. 9B, compared 
with what we found with the K-matrix fit. The fact that both determinations overlap is a good test of 
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Figure 9B. Comparison of the fits to ryg"': polynomial [Eq. (B.3)] given by the dotted lines, and 
with the K-matrix [Eq. (2.4)] (solid line and dark area). 



the correctness of our assumption, for the K-matrix fit, that the contribution of multiparticle channels is 
comparable to the error of the fit itself. 

The fulfillment of dispersion relations with these polynomial fits is just as good as with the K- 
matrix fit; however, the errors of the K-matrix fit are smaller than what we find with the polynomial fit: the 
fulfillment of said dispersion relations may therefore be considered to be marginally better with the K-matrix 
formalism, which is why we only gave results for the dispersion relations with the K-matrix fit. 
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